## Replicate Figure A2
dat2017$unauthorised <- ifelse(dat2017$auth_rev == "1", 1, 0)
dat2017$noprotest <- ifelse(dat2017$auth_rev == "2", 1, 0)

logit1 <- glm(nonans ~ unauthorised +noprotest + age + gender +education +moscowstp + betteroff, data = dat2017, family = "binomial")

fig_A1B<-dwplot(list(logit1),
                vline = geom_vline(xintercept = 0, colour = "grey60", linetype = 2)) %>% # plot line at zero _behind_ coefs
  relabel_predictors(c(unauthorised = "Unauthorized",  
                       noprotest = "No protests",
                       age = "Age", 
                       genderfemale = "Female", 
                       education = "Education", 
                       moscowstp = "Moscow/StP", 
                       betteroff= "Income"))+
  theme_bw() + xlab("Coefficient Estimate") + ylab("") +
  geom_vline(xintercept = 0, colour = "grey60", linetype = 2) +
  ggtitle("Non-Response, 2017 Survey") +
  theme(plot.title = element_text(face="bold"),
        legend.position = c(0.007, 0.01),
        legend.justification = c(0, 0), 
        legend.background = element_rect(colour="grey80"),
        legend.title = element_blank()) 
fig_A1B

ggsave(height=5,
       width=7,
       "Figures/fig_A1B.png")
